Infinite-randomness critical point in the itinerant quantum antiferromagnet 
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We study the quantum phase transition in the three-dimensional disordered itinerant antiferro- 
magnet by Monte- Carlo simulations of the order-parameter field theory. We find strong evidence 
for the transition being controlled by an infinite-randomness fixed point: The dynamical scaling is 
activated, i.e., the logarithm of the energy scales like a power of the length, implying a dynamical 
exponent of infinity. The probability distribution of the energy gaps is very broad and becomes 
broader with increasing system size, even on a logarithmic scale. 
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Phase transitions in systems with quenched disorder 
are an important topic in statistical physics. The critical 
behavior of systems with quenched disorder can be di- 
vided into three classes, according to the behavior of the 
disorder under coarse graining. In the first class, disorder 
decreases under coarse graining, and the system becomes 
asymptotically homogeneous at large length scales. Tech- 
nically, this means the disorder is renormalization group 
irrelevant, and the transition is controlled by a "clean" 
fixed point (FP). According to the Harris criterion |l| 
this happens if the clean FP fulfills v > 2/d, where v 
is the correlation length critical exponent and d is the 
spatial dimensionality In this first class the macroscopic 
observables are self-averaging at the critical point, i.e., 
the relative width of their probability distributions goes 
to zero in the thermodynamic limit [| . 

In the second class, the system remains inhomogeneous 
at all length scales with the relative strength of the in- 
homogeneities approaching a finite value for large length 
scales. These transitions are controlled by renormaliza- 
tion group FPs with finite disorder. Macroscopic observ- 
ables are not self-averaging, the relative width of their 
probability distributions approaches a size-independent 
constant |^|, |) . Examples of critical points in the second 
class include the dilute three-dimensional Ising model, 
classical spin glasses, and various other thermal critical 
points in disordered systems. 

The third possibility occurs when the relative mag- 
nitude of the inhomogcncitics increases without limit 
under coarse graining. The corresponding renormaliza- 
tion group FPs are called infinite-randomness FPs. At 
these FPs the probability distributions of macroscopic 
variables become very broad (on a logarithmic scale) 
with the width increasing with system size. Conse- 
quently averages will be often dominated by rare events, 
e.g., spatial regions with atypical disorder configura- 
tions. Infinite-randomness critical points have mainly 
been found for quantum phase transitions since the disor- 
der, being perfectly correlated in (imaginary) time direc- 
tion, has a stronger effect for quantum phase transitions 
than for thermal ones. Examples include the random 
transverse-field Ising and Potts chains ||, |(| and the 



two-dimensional transverse- field Ising model J?], |§| . 

A natural question is, how general is the occurrence of 
infinite-randomness FPs in disordered quantum systems. 
One prototypical and particularly controversial transi- 
tion is the antiferromagnetic quantum phase transition of 
disordered itinerant electrons in three dimensions which 
is believed to be important, e.g., for a variety of heavy- 
fermion materials H . It has been investigated by various 
methods but no definite results could be achieved. In 
the absence of disorder this transition is controlled by 
a Gaussian FP with mean-field static critical exponents 
and a dynamical exponent of z = 2 H]. According to 
the Harris criterion [Q this clean FP is unstable with re- 
spect to disorder, and the transition, if any, must be in 
class two or three of the above classification. Within the 
conventional perturbative renormalization group [p"lf one 
finds a finite-disorder FP suggesting that the transition 
belongs to class two. However, the renormalization group 
flow diagram, Fig. |l|, shows that a system with weak ini- 
tial (bare) disorder is taken to large disorder at interme- 
diate stages of the renormalization before spiraling into 
the FP. This casts serious doubts on the validity of the 
perturbative approach. Indeed, by taking into account 
the effects of rare regions it was later shown [[l2) that 
the conventional FP is unstable and the renormalization 
group flow is towards large disorder in all of the physi- 
cal parameter space. The ultimate fate of the transition, 
however, could not be determined within this approach. 
Possible scenarios included a complete destruction of the 
transition, a conventional FP (i.e., class two) inaccessible 
by perturbative methods, or an infinite-randomness FP. 

In this Letter we contribute to the solution of this puz- 
zle by reporting results from large-scale computer simu- 
lations of the antiferromagnetic quantum phase transi- 
tion of disordered itinerant electrons in three dimensions 
with uniaxial (Ising) symmetry. We find strong evidence 
for a sharp phase transition which is controlled by an 
infinite-randomness FP similar to those in the random 
transverse field Ising models 0, ||, [| |?|, ||. Specifically, 
we find that the probability distribution of the order pa- 
rameter susceptibility \ (the inverse energy gap) becomes 
broader with increasing system size, even on a logarith- 
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flow to large disorder 




clean Gaussian FP 



FIG. 1: Schematic of the perturbative renormalization group 
flow. A is the disorder strength and u is a measure of the 
quantum fluctuations 



mic scale. The entire distribution scales with ln% ~ 
with tp w 0.39. This implies a dynamical exponent of 
z = oo. The average gap and the average order parame- 
ter susceptibility follow stretched exponentials with dif- 
ferent powers. In the rest of the paper we outline our 
calculations and discuss the results in detail. 

The simulation of the full disordered interacting elec- 
tron problem is numerically very expensive. Conse- 
quently, the system sizes are so severely restricted that 
an investigation of the critical behavior is essentially im- 
possible in three dimensions. Rather than simulating the 
full problem we have therefore performed simulations of 
the Landau-Ginzburg- Wilson (LGW) theory of the phase 
transition which contains only the long- wave-length, low- 
frequency fluctuations of the order parameter. This ap- 
proach is valid close to the transition as long as the (bare) 
disorder is weak so that local phenomena like the forma- 
tion of localized moments do not play an important role. 

The LGW theory of the antiferromagnetic quantum 
phase transition of itinerant electrons in d dimensions is 
equivalent to a classical ferromagnet in d + 2 dimensions 
with the disorder being uncorrelated in the d "space-like" 
dimensions but completely correlated in the two "time- 
like" dimensions jy], [l| . In the following we concentrate 
on three spatial dimensions and uniaxial symmetry. For 
an efficient simulation we remap the field theory to a lat- 
tice model. The quenched disorder is introduced via site 
dilution in the "space-like" dimensions, i.e., the impu- 
rities arc two-dimensional holes in the five-dimensional 
lattice. The classical model Hamiltonian reads 



H 



E 



S(r,T)S(r',T') 



(1) 



Here r = (x,y,x) and r = (ti,T2) are the "space-like" 
and "time-like" coordinates of lattice sites, respectively; 



and the sum runs over all pairs of nearest neighbors. 
S(r,r) = ±1 is an Ising spin, and e r is a quenched 
random variable with values and 1 and an average of 
[e r ] = p (the symbol [•] denotes the disorder average). All 
the results in this Letter are for p = 0.8 while the perco- 
lation threshold for three-dimensional site percolation is 
p = 0.31. The rather weak disorder allows us to observe 
the crossover from the clean Gaussian FP to the infinite- 
randomness FP using the system size as a renormaliza- 
tion cutoff. In the model Hamiltonian the transition is 
tuned by changing the classical temperature T c i which is 
different from the physical temperature of the quantum 
system; the latter is encoded in the length of the system 
in time-direction. 

For the Monte-Carlo simulations we have employed the 
Wolff cluster algorithm |Q. We have studied systems 
with linear sizes up to L = 29 in space direction and 
also L T = 29 in time direction (the largest system hav- 
ing 20 million sites). We found that no more than 200 
sweeps were required for equilibration. Between 1000 and 
5000 disorder configurations were considered, depending 
on system size with at least 1000 measurement sweeps 
per disorder configuration. 

To get an overview over the behavior of the system we 
have computed the average Binder parameter 
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Here M is the magnetization and (•) denotes the thermo- 
dynamic average for a single sample, [g] has the expected 
finite-size scaling form 

[g] = g{tL 1/v \L T /L Z ) conventional FP, (3) 
[ g ] = g(tL ljv ,\nL T /L' 4 ') infinite-randomness FP (4) 

Here t = (T c j — T c )/T c is the distance from the critical 
point. The dynamical scaling is of power-law type at a 
conventional FP but activated at an infinite-randomness 
FP |Q. As a result of its scale dimension being zero, [g] is 
easily analyze: provided L T is scaled appropriately with 
L, a point in parameter space where [g] is independent of 
the system size corresponds to a renormalization group 
FP. In principle, both v and z (or tp) can be determined 
from the finite-size scaling of [g] . 

In Fig. ^ we show the Binder parameter [g] as a func- 
tion of the classical temperature T c i for different system 
sizes L = L T . For small sizes, L — 7 . . . 10, the curves 
cross at T c \ w 7.88. However, for larger sizes this cross- 
ing point becomes unstable and a new crossing emerges at 
T c i ~ 7.96. A comparison with the flow diagram in Fig. 
[I] suggests that the former crossing point corresponds to 
the clean Gaussian FP. (A system with small bare dis- 
order first approaches the Gaussian FP under renormal- 
ization before it goes to large disorder.) To check this 
hypothesis we have performed a series of calculations at 
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classical temperature T d classical temperature T, 



FIG. 2: Average Binder parameter [g] as a function of classical 
temperature T c i for different system sizes. 



FIG. 3: Log. plot of the Binder parameter [g] close to the 
critical point. At T = 7.96 the standard deviation is about 
one symbol size. Inset: Finite-size scaling plot using v = 1. 



T c i = 7.88, independently varying L and L T . The max- 
imum of [g] as a function of L T occurs at L « L T as 
expected at the Gaussian FP. A standard finite-size scal- 
ing analysis based on (j^) yields exponents of v s» 0.55 
and z « 1.0. Since the range of available sizes at this 
unstable crossing is very limited the agreement with the 
expected values vq = 0.5 and zq = 1 is surprisingly good 
jL5|. We conclude that the crossing point at T c ; « 7.88 
indeed corresponds to the clean Gaussian FP. 

We now turn to the crossing point at T c \ as 7.96 which 
we attribute to the critical FP of the phase transition. 
We have again performed a series of calculations inde- 
pendently varying L and L r . However, the dependence 
of [g] on L T for fixed L turned out to be very weak. 
(At T c i = 7.96 and L — 18, [g] remained constant at 
g c = 0.048 within our statistical error of about 6% for all 
L T between 12 and 29.) This can be explained by the fact 
that the critical value, g c = 0.048, is very small, which 
corresponds to weak average correlations. This tends to 
reduce the influence of the boundary conditions and thus 
the sample shape [|l6]. Consequently, finite-size scaling 
of [g] is not an efficient method to study the dynamical 
scaling of this critical point. For all further calculations 
we have therefore used samples with L = L T . Fig. || 
shows the average Binder parameter [g] for systems with 
L > 12 in the vicinity of the critical point at T c i — 7.96. 
The inset shows that the data scale reasonably well with 
a correlation length exponent of v = 1.0. 

In order to further analyze the critical behavior we now 
consider the probability distribution of the order param- 
eter susceptibility Xi i- e -i the inverse energy gap. At a 
conventional finite-disorder FP the distribution of x/[x] 
should be size-independent (|, ||]. In Fig. |^ we show 
this distribution at the critical point, T c i — 7.96, for 
different system sizes. Clearly the distribution becomes 
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FIG. 4: Distribution of the order parameter susceptibility 
x/[x] a t th e critical point T c i = 7.96 for different sizes (at least 
2300 disorder configurations were used for the distributions). 
Inset: Dependence of A^/[x] on L. 

broader with increasing L. This can also be seen from 
the inset which shows the relative widths of the distribu- 
tion as a function of system size. This suggests that the 
phase transition is controlled by an infinite-randomness 
FP rather than a conventional one. 

At an infinite-randomness FP we expect activated scal- 
ing 0], lnx ~ L^. In Fig. [| we show a corresponding 
scaling plot of the distribution of ln^ The distribution 
scales with an exponent of ip = 0.39. The inset shows the 
size-dependence of the logarithm of the typical suscepti- 
bility, In Xtyp — [ m x]- The solid line is a power-law fit, 
In 

Xtyp ~ L^, which gives the same exponent of ip = 0.39. 
We now consider the average susceptibility [x] and the 
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FIG. 5: Distribution of x = hix/L* at T ci = 7.96 using 
V> = 0.39. Inset: Dependence of \n\t yv on L. The solid line is 
a power law fit with the same exponent ip — 0.39. 



an infinite-randomness FP for the XY and Heisenberg 
cases, too. However, in Ref. || it was found that the 
random-singlet FP in localized XY or Heisenberg anti- 
ferromagnets |l7], [fj|, [l9| , another example of an infinite- 
randomness FP, becomes unstable for d > 1, suggesting 
that a continuous order parameter symmetry tends to 
weaken infinite-randomness FPs. It thus remains unclear 
whether the antiferromagnetic quantum phase transition 
of disordered itinerant electrons with XY or Heisenberg 
symmetry is controlled by any of the know types of FPs 
or whether it belongs to a new class. 
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CLIC parallel supercomputer of Chemnitz University of 
Technology, using a total of about 6000 CPU days. 



average gap [1/x]- Since for large L the distribution of 
lnx becomes very broad, [x] will be dominated by the 
large-(lnx) tail of the distribution. For large x = h\\l 
the distribution in Fig. || approximately falls off as e~ c * x . 
For large L the integral for [x] can be treated in sad- 
dle point approximation giving an asymptotic size de- 
pendence of ln[x] ~ L 2 ^ . The average order parameter 
susceptibility increases faster with system size than the 
typical one. Analogously, the average energy gap [1/x] 
will be dominated by the small- (In x) tail of the distribu- 
tion. However, the accuracy of our data is not sufficient 
to determine the functional form of this tail. For system 
sizes L = 12 to 29 the susceptibility distribution is not 
sufficiently broad yet for the above saddle-point argu- 
ment to apply. In this size range we numerically obtain 
the effective relations ln[%] ~ L QA5 and — ln[l/x] <~ L 37 . 
While the exponent of ln[%] is larger than if), it is much 
smaller than the asymptotic value of 2ip = 0.78. 

In summary, we have found strong evidence for an 
infinite-randomness FP at the antiferromagnetic quan- 
tum phase transition of disordered itinerant electrons in 
three dimensions. This FP has properties similar to those 
found in the random transverse-field Ising model in one 
and two dimensions. We conclude with two remarks. 
First, the systems we were able to simulate are rather 
modest in linear size because the effective dimensionality 
of the LGW theory is very high. Therefore all quantita- 
tive results for critical exponents should be understood 
as effective rather than asymptotic values. Second, the 
results obtained here are for the case of uniaxial sym- 
metry. The qualitative features of the perturbative flow 
diagram in Fig. [I] including the runaway flow to large 
disorder, are the same for Ising, XY or Heisenberg sym- 
metry. This suggests that the transition is controlled by 
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